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Abstract 

The slow-to-start mechanism is known to play an important role in the particular 
shape of the Fundamental diagram of traffic and to be associated to hysteresis ef- 
fects of traffic flow. We study this question in the context of exclusion and queueing 
processes, by including an asymmetry between deceleration and acceleration in the for- 
mulation of these processes. For exclusions processes, this corresponds to a multi-class 
process with transition asymmetry between different speed levels, while for queueing 
processes we consider non-reversible stochastic dependency of the service rate w.r.t the 
number of clients. The relationship between these 2 families of models is analyzed on 
the ring geometry, along with their steady state properties. Spatial condensation phe- 
nomena and metastability is observed, depending on the level of the aforementioned 
asymmetry. In addition we provide a large deviation formulation of the fundamental 
diagram (FD) which includes the level of fluctuations, in the canonical ensemble when 
the stationary state is expressed as a product form of such generalized queues. 

1 Introduction 

In the study of models for traffic [1], properties of the FD, which gives a relation 
between the traffic flux and the density or alternatively the dependence between the 
speed and the flux, or the speed and the density, play an important role. In the three 
phases traffic theory of Kerner [2] , the traffic phase diagram on highways consists of 
the free flow, the synchronized flow and the wide moving jam. In the free flow regime, 
at low density, the flux is simply proportional to the density of cars; in the congested 
one, at large density, massive clusters of cars are present, and the flux decreases more 
or less linearly with this density; in the intermediate regime, the relation between flux 
and density is largely of stochastic nature, due to the presence of a large amount of 
small clusters of cars propagating at various random speeds. It is not clear however 
whether in this picture these phases, and especially the synchronized flow phase, are 
genuine dynamical or thermodynamical phases, meaningful in some large size limit in 
the stationary regime, or are intricate transient features of a slowly relaxing system. 
In fact there is still controversy about the reality of the synchronized phase of Kerner 
at the moment [3]. 

In order to analyze properties of the FD from a statistical physics perspective, 
we look for a simple extension of the totally asymmetric exclusion process (TASEP) 
as well as the zero range process able to encode the fact that vehicles in the traffic 
may accelerate or brake. Some asymmetry between these two actions is empirically 
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known to be responsible of the way spontaneous congestion occurs. This is referred as 
the slow-to-start mechanism, not present in the original cellular automaton of Nagel- 
Schreckenberg [3], but in its refined version like the velocity dependant randomized 
one (VDR) [5] which exhibits a first order phase transition between the fluid and the 
congestion phase and hysteresis phenomena [6j associated to metastable states. So a 
realistic stochastic model should have in particular the property that some spontaneous 
symmetry breaking among identical vehicles can occur, as can be seen experimentally 
on a ring geometry for example [7]. With such a model at hand we would like to 
provide a method to compute the FD and study emergence of non-trivial collective 
behaviors at macroscopic level. 

The paper is organized as follows: in Section [2] we start by defining a simple and 
somewhat minimal generalization of TASEP and describe some of its basic properties. 
In Section [3] a family of zero-range processes is introduced which service rate depends 
stochastically on the state of the queue and on which there are relevant mappings of 
multi-speed TASEP; we determine sufficient conditions for having a product form for 
the invariant measure of such processes coupled in tandem. In Section [4] we provide a 
large deviation formulation giving the FD along with fluctuations on the ring geometry 
in the canonical ensemble, when the steady state has a product form. Finally, in 
Section[5]we analyze the jam structure of a generalized queueing process corresponding 
to the model introduced in Section [2] and interpret direct simulations of the multi- 
speed TASEP in this light. 

2 Multi-speed exclusion processes 
2.1 Model definition 

The model we investigate is a multi-type exclusion process, generalizing the simple 
exclusion process on the line, introduced by Spitzer in the 70's [H1[H] 5 combined with 
some feature of the Nagel-Schreckenberg cellular automaton [4]. 

In the totally asymmetric version of the exclusion process 8, 9 (TASEP), particles 
move randomly on a 1-d lattice, always in the same direction, hopping from one site 
to the next one within a time interval following a Poisson distribution and condition- 
ally that the next site is vacant. In the Nagel-Schreckenberg cellular automaton, the 
dynamics is in parallel: all particles update their positions at fixed time intervals, but 
their speeds are encoded in the number of steps that they can take. This speed can 
adapt stochastically, depending on the available space in front of the particle. 

The model that we propose combines the braking and accelerating feature of the 
Nagel-Schreckenberg models, with the locality of the simple ASEP model, in which 
only two consecutive sites do interact at a given time. The trick is to allow each car 
to change stochastically its hopping rate, depending on the state of the next site. For 
a 2-speed model, let A denotes sites occupied by fast vehicles, B sites occupied by 
slow ones and let O denote empty sites; the model is defined by the following set of 
reactions, involving pairs of neighbouring sites: 



AO 




■OA 


simple move of fast vehicle 


(1) 


BO 


f ''\ 


OB 


simple move of slow vehicle 


(2) 


BO 




■AO 


slow vehicle spontaneously accelerates 


(3) 


AO 


6 


■ BO 


A brakes behind O = A or B 


(4) 



fi a , Hb, 7 and S denote the transition rates, each transition corresponding to a Poisson 
event. The dynamics is therefore purely sequential as opposed to the dynamics of the 
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Nagel-Schreckenberg model. It encode the fact that a slow vehicle tends to accelerate 
when there is space ahead, while in the opposite case, corresponding to Q it tends to 
slow down. The main mechanism behind congestion, namely the asymmetry between 
braking and acceleration is potentially present in the model when 7 is different from 
S. Our model is in fact similar to the model of Appert-Santen [10], in which there is a 
single speed, but particles have 2 states: at rest and moving with possible transitions 
between these 2 states. 

To define fully the model, its boundary conditions have to be specified, either 
periodic in the ring geometry or with edges. In the latter case, additional incoming 
and outgoing rates have to be specified, depending on whether we want to model a 
traffic light a stop or simply a segment of highway for example. In this paper we will 
consider only the ring geometry. In the following we will refer to this specific model 
as the acceleration-braking totally asymmetric exclusion process (ABTASEP). 



2.2 Known special cases 

Let us first notice that this model contains and generalizes several sub-models which 
are known to be integrable with particular rates. The hopping part 1 1|2 1 of the models 



is just the TASEP when fj, a = fit, which is known to be integrable. Its generaliza- 
tion to include multiparticle dynamics with overtaking is the so-called Karimipour 
model [111 I12j which turns out to be integrable as well. Matrix Ansatz method al- 
lows in some cases, like the ASEP with open boundary conditions, to describe the 
stationary regime of the models, using the representations of the so-called diffusion 
algebras [13]. The acceleration/deceleration dynamics is equivalent to the coagula- 
tion/decoagulation models, which are known to be solvable by the empty interval 
method and by free fermions for particular sets of rates [14] . but the whole process is 
presumably not integrable. 



2.3 Relation to tandem queues 

In some cases, the model can be exactly reformulated in terms of generalized queueing 
processes (or zero range processes in the statistical physics parlance), where service 
rates of each queue follows as well a stochastic dynamics [TS] . In this previous work we 
however considered exclusion processes involving three consecutive sites interactions 
in order to maintain the homogeneity of labels in particle clusters. The mapping works 
only on the ring geometry, by identifying queues either with 

(i) cars: clients are the empty sites. 

(ii) empty sites: clients are the vehicles; 

In our case the mapping of type (i) is exact. In the corresponding queueing process, 
queues are associated either with fast or slow cars, having then service rates jj, a or /it,. 
Slow queues become fast at rate 7 conditionally to having at least one client, while 
empty fast queues become slow at rate S. 

The mapping of type (ii), is more informative with respect to jam distribution, 
but is not possible with transitions limited to 2-consecutive sites interactions, because 
in that case homogeneity is not maintained in the clusters, and additional information 
to the number of cars and the rate of the car leaving the queue is needed to know the 
service rate of the queue. For this mapping we have to resort to some approximate 
procedure as shall be exemplified in Section [5] 

Another way to circumvent this would be to start from a a slightly different defi- 
nition of our initial model, by not attaching speed labels to cars, but instead to empty 
sites. Let A and B denote empty site with respectively fast and slow speed, and V 
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denote sites occupied with a vehicle, the set of transitions reads then 



'VA 




AV 


simple move on fast site 


VB 


l-> b 


BV 


simple move on slow site 


VA 


6 


VB 


fast site become slow 


VB 


1 


VA 


slow site become fast 



(5) 



where V is an unoccupied site i.e. either A or B. The mechanism in this model is 
that the empty sites visited more recently have a slower associated speed than others, 
leading possibly to congestion instability. The mapping to tandem queues is then 
straightforward: queues are associated to empty site of type A or B with corresponding 
service rate fi a or pb\ fast queues with at least one client have a probability per unit 
of time S to become a slow queue, while slow queue which are empty become fast with 
probability rate 7. This model generalizes directly to any speed levels. Interestingly, 
one see that in such a model the service rate of queues are somehow related to the 
lifetime of clusters. In this paper we will focus on model of Section [2. 1| deferring for 
future work the study of model Q . 

2.4 Numerical observations 



Based on numerical simulations on the ring geometry, we make some observations 
concerning the phenomenology of the model, depending on the parameters. This is 
illustrated on Figure [I] When no asymmetry between braking and accelerating is 




Fig. 1: Space-time diagram with 2 speed levels (a)-(c) and 3 in (d). Time is 
going downward and particles to the right. Red, green and blue represent 
different speeds in increasing order. The size of the system is 3000 except 
for (b) where it is 100000. Setting are p a = 100, p b = 10, 7 = 10, 6 = 1 
for (a) and (b) and S — 10 for (c), all with density p = 0.2. 

present, as in TASEP on a ring, no spontaneous large jam structure is observed. As 
the density p of cars increases, one observes a smooth transition between a TASEP 
of fast particles for small p to a TASEP of slow particles around p ~ 1. Instead, 
when the ratio Sf-y is reduced, there is a proliferation of small jams. Below some 
threshold of this ratio, we observe a condensation phenomenon, associated to some 
critical value of the density: above this critical value, after a time which presumably 
scales as some power of L, one or more large jams, which absorb a finite fraction 
of the total number of cars may develop. If there are many of such large jams, a 
competition occurs, which end up possibly in the long term with one single large 
condensate depending on the relative value of 7 and pt and also of the size N of the 
system. It is tempting to interpret this as a condensation mechanism at equilibrium 
in the canonical ensemble |16j . by combining the Nagel-Paczuski [17] interpretation 
of competing queues with some results in a previous work rT5] , which in the context 
of tandem queues on a ring, allows this condensation mechanism to take place if the 
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apparition of slow vehicle is a sufficiently rare event. As such, the mechanism for 
congestion is then understood qualitatively as follows: assume that a jam of size s 
appears, created somewhere by a mutation A — > B. Let t(s) be the expected waiting 
time in this queue. The probability that a vehicle is still of type A when it reaches 
the end of the queue is simply 

p(n = Ha) = e~ st{s) = p s . 

Assuming that s does not change much during t(s), then we have 

t{ S ) = s(—p s + M -Ps)), 

if we take into account the fact that the vehicle may accelerate at rate 7 before leaving 
the queue. Anyway, this gives t(s) self-consistently and yields qualitatively that, at the 
beginning of the process, there is a distribution of jams with various effective service 
rates. As time evolves, either long-lived jams with effective service rates fJ. a {/J-b + 
7)/(/i a + 7) are able to survive, then new small jams can never appear because their 
effective service rates are strictly larger, so in the end there is a competition between 
the existing jams; the largest one eventually remains alone after waiting a large amount 
of time. Typically this happens when 5 <C /it and 8 < 7 (Figure |TJb) ) . If the situation 
with a single jam is not stable, then no large jam may develop at all, and only small 
fluctuations are to be observed. When S > 7 one likely observes the kind of jams seen 
in Figure [TJc). 



N=10000 va=100vb = gamma=delta=10 N=10000 va^lOO vb=gamma=10 delta=l 




Traffic density Traffic density 



Fig. 2: Fundamental diagram with 7 — 5 (left) and 7 = 105 (right) with small 
spontaneous rate of emission and escape of particles. Trajectories of 
$1 = X^q 1 VaAiOi+i + HbBiO i+ i and <J> 2 = Y^Jq VaM + VbBi are 
displayed. 

The effect of the asymmetry is also clearly observed on Figure [2] where we allow 
particle to enter or quit the system with some very low probability rate compared 
to others. The global density of cars then performs a random walk, and by looking 
at trajectories in the FD plane, we see hysteresis effect for 5 < 7. We have also 
simulated a model with 3 speed levels, see Figure [TJd). In that case, small jams may 
have different speeds, depending on which type of slow car is leading. Then a cascade 
mechanism takes place, slow speed regions generate even slower speed clusters of cars 
and so on, and some kind of synchronized flow is observed. 

To conclude this section, let us finally remark that Figure [TJa) is very reminiscent 
of coagulation-decoagulation process, which is somewhat expected from the previous 
discussion, since it is present in the equations. One could therefore expect the jam 
structure to share some properties with the directed percolation process taking place 
on the ring geometry. 
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2.5 Hydrodynamic equations and solitons 

Although it is not clear whether the hydrodynamic limits is valid in our context, one 
can at least write down the corresponding equations on the density p a (x, t) and pb{x, t), 
with x = i/N and N large, expected in such a limit: 

dtpa + V a d x (p a p~) = -5p a p + 7 ' pbp (6) 

d t pb + Vbd x (pbp) = Sp a p — 7/0i,p (7) 
with p — p a + Pb the total density of cars, p — 1 — p the density of empty sites and 

(JV) 

dot .. M a ,f> 

v a ,b = lim 



AT->oo N 

if the rates p^ are allowed to be rescaled when N is varied. This is an hyperbolic 
system of equations since the matrix 



M 



has two real eigenvalues 



V a (l-2p a — pb) —V a pa 

— V b pb V b (l — p a — 2p b ) 



A* = ^ (v a + v b - v a (2p a + p b ) - v b (p a + 2p b ) ± V~A, 



with 



A = (u (l - 2p a - pb) - Vb(l - p a - 2pb)) +4v a v b p a p b > 0. 

The method of characteristic applied to this system is however not giving much insight 
into it, because of the non-linear shape of the characteristics. At steady state we can 



P = i 



I'nr/n 




Fig. 3: 

nevertheless look for backward travelling waves solutions, of the form p a , b (x + vt), 
with v the propagation velocity. Replacing dt by vd x in ( 6j6 ( and summing up the 
two equations one can check that the flux 



is a conserved quantity 



<J> = (V a pa + V b pb)p + Vp, 



dx 



= 0. 



This allows to express p a and p b as function of p: 
<J> — pv — Vbpp 



Pa = 



and 



(v a ^ Vb)p 6 (Vb - v a )p 

After rescaling w.r.t v a the model has only 3 independent parameters: 



$ — pv — v a pp 



def Vb dot S 

V = — ) K = ~, 

V a 7 



dof 7 
v = — . 

Vb 



(8) 



(9) 
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and we rescale as well the soliton speed v = f v/v a and the flux <j> == &/v a - Replacing 
p a and pb with their expressions in J9J, we obtain an equation relating 9a; p and p: 

da; P 0(p)' 

with 

C=£(1-„k) 
and where P and Q are both polynomials of degree 3: 

P(p) =' p(p - p- ) (p - p + ) + C- (1 - K)p + KV) (<j> - V) 

Q(p) = p 2 (p - pi) + (0 - ? 



where 



2rj 



,1 + 7? ± 1 + (1 - k)v ± \/A 

Pi = H v, p - 



r, ' 2(1 -t]k) 

are the roots when <f> = v with the discriminant 

A = (1 - (1 + k)v) 2 + 4(77 - v)kv, 

assumed to be positive. Assuming n < 1 we have p _ £ [0, 1] and p + > 1 [resp. 
p + £]p~, 1]] when v < rj [resp. w > 7;]. 1/C gives the length scale of the interface 
between regions of different densities. In particular if p a ,b are maintained at a fixed 
value the interface becomes a step function when N is large. 

Let pmax and p m ; n the extremal values of the density. We expect a solution which 
decreases from p max to p min with R(p) staying strictly positive for p G [pmin , pmax] ■ 
Since R(p) is a single valued function of p, there must be a discontinuity, i.e. a shock 
to relate these two values as depicted on Figure [3] We don't expect everything to be 
discontinuous at this point, in particular, as long as the shock is not at rest (v = 0) the 
proportion of fast cars u = f p a /p should be continuous, since cars are slowing down at 
finite pace 5. This means that pmin and p m ax are not independent. Instead, letting 
u* the value of u at the shock, from j9||, they are both solution of (u* == 1 — u*) 

(v a u* + v b u*)p 2 — (v a u* + v b u* + v)p + = 0. 

Given that the sum and product of the two roots have to be smaller than 2 and 1 
respectively, we deduce that v and <f> are both strictly smaller than v a u* + VhU*. In 
addition the constraint p max < 1 implies <f> > v. As a consequence, for p g [0, 1[ we 
have Q(p) > 0, vanishing only at p — 1 in the limit case $ = v. Concerning P(p) it has 
then 3 roots pi,i £ {1, 2, 3} with < pi < P2 < 1 < pz- To obtain a general solution, 
there are two free parameters, v and pmin and therefore 2 constraints, associated to 
fixing the length of the system to one the mean density to po. These read 



/' 



d ?Q^=C and r X dp^m = C P o. 



P P (P) J Pmi „ P p (p) 



When condensation occurs we expect to observe a full congestion (p = 1) at some 
point. This entails from (JsJ> that cj> = v, yielding in that case the equation: 

d P _ C (P + - P)(P~ P~) 



dx p — pi 
and 

u - VP r , -, 

pa = for p < 1. 

1-77 

On one hand, we should have p + > 1 to be able to have a point where p = 1, 
consequently v is necessarily smaller than 77; on the other hand p a becomes negative 
when p — ► 1 if v < 77. As a result, there is no consistent way to obtain genuine 
condensation in this hydrodynamic framework. 
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3 Tandem queues with stochastic service rates 

3.1 Definition of the process 

As already mentioned in Section |2.3| in some cases, on the ring geometry, exclusion 
processes can be mapped onto a queueing process. These have a fixed number of 
queues, organized in tandem and in each queue, both the number of clients and the 
service rate are stochastic, the service rate being not necessarily a function of the 
number clients. In some case, each queue taken in isolation may not be a reversible 
process, with some hysteresis feature relevant to traffic. This family of models is 
defined more formally as follows: let Q — (A/ - , C) a network of queues with dynamical 
(stochastic) service rates. By dynamical service rates, we mean that each single queue 
i £ Af is represented by a vector Zj(t) = (nj(£), //»(£)) £ E t C N x R + . n,i(t) is the 
number of clients and n%{t) is a service rate, which represents the global transition 
rate from Zi to z[ = (rii — l,/4) £ (zi) (V~(z) is the set of points in Ei having 
one client less than z). This process is a generalization of queueing processes with 
stochastic service [181 1151 [14] . From the viewpoint of the statistical physics it is an 
extension of the zero-range process, obtained by adding internal dynamics. 

Two sets of transition probability matrices p^(z, z') and one set of transition rates 
qf(z,z') are introduced to complete the definition of the process. When a client get 
served in queue i, there is a transition from the state of the departure queue Zi to 
one of the state z' £ V~ (z) with a rate [ii(t)p~ (z, z'); in the destination queue, a 
transition occurs from the state Zi+i into one of state z' £ V i '^ 1 (z) with probability 
pf +1 (z, z'). We have the normalizations, 

E Pi (*,*/) = !■ (10) 

Additional internal transitions are allowed, where the service rate fn of queue i changes 
independently of any arrival or departure. The intensities of these transitions are given 
by the set q®(z,z'), z' £ V°(z) of transition rates (V^z) is the set of points in Ei 
having the same number of clients as z). The combined set of transition rates, 

qi{z,z) = \pf(z,z')l {z , €V + iz)} + n(z)pi (z, z')l {z , ev - (z)} +qi{z,z')l {z , ev o {z)) , 

defines for each i £ J\f a continuous time Markov process representing the dynamics 
of one queue taken in isolation with arrival rate A. The joint process is then entirely 
specified by the state graph of the single queue which set of nodes {(n, fj,)} is a subset 
of N x R+, and in which each transition is represented by an oriented edge (see e.g. 
Figures |i]and[5|. 

3.2 Product form of clusters at steady-state 

In the stationary regime, the invariant measure of a network of reversible queues is 
known to factorize into a product form in general [19] . which allows to compute many 
equilibrium quantities explicitly. With a dynamical rate it is interesting to check for 
such properties and in this section we study the conditions under which the stationary 
state has a product form. For simplicity we restrict the analysis to the case where Q 
corresponds to a tandem of N queues where the clients leave queue i to enter in queue 
(i + 1) mod (N), for any i = 0, . . . , N — 1. We can establish the following sufficient 
conditions: 

Theorem 3.1. Let denotes the steady state probability corresponding to queue i 
taken in isolation and fed with a Poisson process with rate X. If the following partial 
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balance equations are satisfied (see FigureUty, 

E V( Z )P7( z > z i) n i( Z ) =^i( Z i)> 
fl(zi)lTi( z i) + 9i( Z iT Z )^H Z i) = 



(11) 



(12) 



i/ien t/ie joint probability measure of the network has the following product form at 
steady state: 



P(S = {zi,i<EAf}) = 



Proof. Using ||10|), the global balance equations reads 



E 



(13) 



/ l\ —I I \ + I II \ ^ii 2 ) 7v i+l{ z ) 

y. v( z >Pi ( z , z i)Pi+i( z \j \ x / — ^ 



+ E 



*H Z ) 



P(S = {zi}) 



= E[^)+ E q°(*i,*j\P(S = {zi}) 

ieM z£V°(zi) 

To find a sufficient condition select term i in the first term of the left hand side and 
terms i + 1 otherwise in the preceding equation, yielding Vi £ A/", 

E v( z ')Pi( z ',^)pt+i( z ' , , z i+i)'Xi( z ') 7r i+i( z '')+ 

^ <? l °+i(2,2: l+ i)7r l A (2i)7r l A +1 (2) 
^evP +1 (z i+ i) 

= pi(z l+ l)7r i A (2 l )7r i A +1 (z i + l) + ^ g i °+l(2: l +l,z)7I- l A (z0 7r l A +l( 2: i + l)- 

After multiplying this last equation by A, this sufficient condition rewrites 

E ^{ Z ')Pi{ Z ' -: Z i)^i{ Z )) ( A E Pt+l( Z '' Zi + l)^+l( 2 ")) 



z"ev i+1 ( z , + i) 



= A7Tj (Zi)X 

^(2 i+ i)7r A +1 (2; i+ i) + ^ g l +1 (z l+ i,z)7r A +1 (2: l+1 ) - (z, 2 4+ i)tt A +1 (z)J , 
»£Vi C +1 (»i+i) 



which is satisfied Vi £ A/" under the two sets of sufficient conditions (111 and ( 12 1. ■ 

Note that reversible processes are special cases of processes obeying ( 1 1| 12 l and, 
in this respect, our result is an adaptation to our context of the general results of 
Kelly concerning product forms in queueing networks [TS] . The next obvious question 
is then whether there exists any such process which is non-reversible. 
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3.3 Examples 

Let us give two examples of non-reversible queues yielding a product form. The 
isolated-queue state-graph corresponding to the first one is sketched in Figure |4jb). 
Using the labelling (0) = (0,Ho), (1) = (l,A*i), (2) = (l,/i 2 ), (3) = (2, M4 ) and 
(i) = (n — 2, in) for i > 4 for the different states of the isolated queue, these partial 
balance conditions written for the nodes (i),i < 4 yield four independent relations out 
of six, between the single queues steady states probabilities 71^: 

XlVo = jUl7Tl + ^27T2, A7Ti = ^3P3l7T3, 

g2l7T2 = (<7l2 + Ml) 75 "!! ^2 = ^3P3^3, 

which are actually compatible iff: 

P32 _ gX2 + Ml 
P31 ?21 

Letting p = P31 = 1 — P32, this corresponds to 

921 

P = • 

Mi + 912 + q2i 

For i > 4, detail balance holds anyway. When this hold the joint measure of the 



tandem queue takes the product form ( 13 1 whereas this example is by construction 



never reversible (transitions — > 1 is absent), with the single queue invariant measure 



3 Tandem queues with stochastic service rates 



11 




simply reading 



7Tl 



Xqi 



Hl(l2 + <?12M2 + 921 Ml 



-7T0, 



712 



A(mi + g 2 i) 



Ml^2 + gi2A*2 + <?2lMl 



-TO, 



T3 



A 2 (/ii + glj + g2l) 

+ gi2A*2 + 921/ii 



-710, 



71n>4 



n 



Aqi2 



MlM2 + 912/12 + 921^1 



-TO, 



where 7io is set upon normalization. 

Our second example is depicted in Figure [5] It represents an exclusion process on 
a ring, where particles, to hop to the next site have to follow a cascading process. For 
sake of simplicity let us limit ourselves to a two stage cascade, where empty site are of 
two different types: either blocking type B and non-blocking type A. The transition 
rules then read: 



V A — i> BV vehicle hopping 

VB V A blocking site becomes non-blocking 

AB AA blocking site becomes non-blocking 



(14) 
(15) 
(16) 



This model is mapped on a composite queueing process as follows: cars represent 
composite queues, with client arriving in the type B (first stage of the composite 
queue) and able to jump to the second stage in type A when they reach the server of 
the queue, i.e. when the left neighbor is of type A. The whole tandem of queues has a 
product form since each single queue taken in isolation is reversible. Considering then 
the queue depicted on Figure|5|b) as an effective equivalent to the composite one, it is 
clearly non-reversible. Nevertheless the partial balance sufficient condition for having 
a product form are fulfilled. Indeed in the example of the figure the isolated effective 
queue has for invariant measure 



7r(n, ft 



Mi + M2 



) = (l-pi)(l- P 2K, 



7r(n, M = Mi) = (1 - Pi)(l - (to) 



Pi - P2 
Pi - P2 



with pi t 2 = A/mi,2- Computing the fluxes yields then the desired partial balance con- 
ditions. It seems that this procedure of replacing composite of reversible queues with 
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effective one yield automatically non-reversible queues fulfilling the partial balance 
conditions. 

The analysis of the fundamental diagram follows then almost straightforwardly 
when such product form is present, although some non-trivial feature are also expected 
depending on the value of the spectral gap between the steady-state and the first 
excited state. 



4 Fundamental diagram for product form distributions 
4.1 Large deviation formulation 

In practice, points plotted in experimental fd studies are obtained by averaging data 
from static loop detectors over a few minutes (see e.g. f5]). This is difficult to compute 
from our queue-based model, for which a space average is much easier to obtain. The 
equivalence between time and space averaging is not an obvious assumption [6], but 
since jams are moving, space and time correlations are combined in some way [17] 
and we consider this assumption to be quite safe. In what follows, we will therefore 
compute the fd by considering the conditional probability measure P{<f>\d) for a closed 
system, where d represents the spatial density of cars and 4> the normalized flux: 



' L number of queues 



N 



N + L' .,, , 
' with < 



N + L' 



L 



N — 'Y_ l rii number of vehicles 

i=i 

L 

$ = Mil{ni>o} integrated flow 



The numbers N of vehicles and L of queues are fixed, meaning in the statistical physics 
parlance, that we are working with the canonical ensemble. If we assume that we are 
in the conditions of having a product form presented in the preceding section, for the 
stationary distribution with individual probabilities 7r A (n, fi) associated to each queue 
taken in isolation then, taking into account the constraints yields the following form 
of the joint probability measure: 

P({m,iH}) = * ( ^zS[ lW<) I>W0, 

i — 1 

with the canonical partition function 

L 

Zi[N] = f Y, •KAT-EirnOn^K.^), 

where S denotes now the usual Dirac function. When <j> is interpreted as a continu- 
ous variable, the properly normalized density-flow conditional probability distribution 
takes the form 

with 

L 

Z L [N,$] d = f *(^-Efcini)*(*-EtiWl{« i >o})n^(n i ,/ii), (18) 

Note (by simple inspection, see e.g. [T5]) that P((f>\d) is independent of A. Zl[N] 
and Zl [N, <3?] represent respectively the probability of having N vehicles and the joint 
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probability for having at the same time N vehicles and a flux <E>, under the uncon- 
strained product form. Under this product form, on general ground, we expect d and 
(j> to satisfy a large deviation principle (see (20] f° r a recent review connecting large 
deviations to statistical physics), i.e. that there exist two rate functions 1(d) and 
J(d, <j>) such that for large L 

Z L {N)^e- LI{d) , 
Z L [N, $] x e - LJ{iM , 

such that we expect a the large deviation version of the fundamental diagram of the 
form 

P(ct>\d) xe- Lif(0|d) . (19) 

with 

K{<j)\d) = J(d,<j>)- 1(d). 

When there is one single constraint like for Zl(N), the large deviation expression can 
be obtained by saddle point techniques like in |16l I21j . For more than one constraint 
it seems easier to work variationally. Consider L n ,j* the number of queues having n 
clients and service rate \i. Assuming fi runs over discrete values, the partition function 
( 18 1 can be recast in: 

z l (n,*)= n L 'l ") £n "' 

x 5(L -J2 L ^) S ( N ~/Z nLn ^) 5 (® ~ zZ » Ln ^) ( 2 °) 

As can be seen explicitly in the following, when there is no condensation, i.e. no 
heavy tail in the L(n, fi) distribution, we are in the condition of the Sanov theorem |2U| 
to get an asymptotic expression of the partition function. Namely, 

Z L (N,*)tz J2 exp(-L^])5(C 1 [ J /])5(C 2 [ J /])5(C 3 M), (21) 

{s/(".m)} 



with 



cf 



L 

the fraction of queues having n clients and service rate fi, J- the large deviation func- 
tional and Ci,2,3[j/] the constraints: 



Ci[y] = ^2y(n,n) - 1, 



Ca[l/] = W(n,v) - y 



.. d 

n>0,fi 



This is obtained explicitly by making use of the Stirling formula in (20 1, keeping only 
the leading terms to get: 

.F(y) = £y(r^)log4^4- (22) 



4 Fundamental diagram for product form distributions 



14 



The rate functions are then obtained using the contraction principle |20) : 
1(d) = inf T[y] 

{»:Cl[»]=0,C 2 [v]=0} 

J(d, 4>) = inf F[y] 

{v-Cily]=0,C 2 ly]=0,C 3 [y]=0} 

which apply here simply because J- is convex and the constraints are linear. The 
variational principle consists in to look for a distribution y minimizing the mutual 
information with the product form while satisfying the constraints; in this sense it is 
equivalent to the standard mean field approximation. A simple generalization would 
then amount to consider, in case where the product form is not valid, the 2-server 
problem to build a joint-law based on pairs 7r (m, 712, H2) of nearest neighbor de- 
pendencies and to use the Bethe approximation. 

Let us introduce the moment and cumulant generating function g and h associated 
to tt a , 

OO 

fl(«j*) = "" A (n,M)e Sn+tM and h(s,t) = log[g{s,t% 

where it is assumed by convention that the rate fi is zero in absence of client. The 
stationary point then reads 

7T>, M ) e " A "+^ 

y{n,n) = — ' , (23) 

with the Lagrange multipliers X n and A M associated to constraints C2 and C3, are 
implicitly given by 

— (Xn(d,<t>), A M (d, tj>)) = 3— j, 

(24) 

J^(A„(d^),A M (d^))= T ^. 

Inserting the variational solution ( 23 1 into J- leads to the resulting expression of the 
rate function as the Legendre transform of h 

J(d, <t>) = j^n{d, 4>) + ^r d x M; 4>) - h(X n (d, 4>), K(d, <t>)) 

which could actually be obtained directly from Cramer's theorem because {(71^,^1;)} 
are iid |20| . Concerning the other rate 7(d) function, if X' n (d) is associated to the first 
constraint, while A M is set to zero we get 

1(d) = j^XUd) - h(X' n (d),0) 

in term of the cumulant generating function h, which completes the determination of 
the large deviation FD 1191). 



4.2 Gaussian fluctuations around the fundamental diagram 

The deterministic part of the fundamental diagram is then obtained for (j>(d) such that, 
dJ{d,4>) \ = A M (d,0(d)) =Q (25) 

d(f> Id, 0(d) 1 — d 

When considering small fluctuations rather than large deviations, we recover a central 
limit theorem version of the fundamental diagram. In that case, replacing in (211 T 
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with its quadratic approximation around some stationary point value y* of y yields 
the following expression for the partitions functions: 



Zl(N) 



exp(-LJ(rf))) 
y/2nLh sa (X' n ,0) 

cxp( — LJ(d, (j>)) 
^(27rL)2det(#*(A n ,A M )) 



(26) 
(27) 



with (A»,A M ) and A'„ solving the constraints (J24J) for a given value of d and and 
denoting 

h ss h s t 
hts htt ' 



(28) 



with the use of shorthand notations for the derivatives, to represent the covariant ma- 
trix between the charges of the queues and the flux. This matrix is actually the Hessian 
in the dual representation, associated to the Legendre transform of the free energy (see 
Appendix |A] for details). We obtain for the large deviation FD the expression 



p(4>\d) = 



L 



27r(l-d) 2 det(ff*(A n ,A„)) 



exp{-LK{4>\d)). 



The Taylor expansion at second order in — 0(d) at fixed density is performed, using 
(25) so that 

J(d,4>(d))=l(d). 



and 



o 2 j(d, <t>) i ay(ri,0) 

d(j)2 UMd) l _ d d(f> UMd) 



h ss {X n , Xfj.) 



d) 2 det(i?*(A„,A^))' 



The last equality is obtained by remarking that from the constraints ( 24 1 we have 

dX n d\fi 
-^r^h ss + ^—hst = 



dX n <9A M 



1-d' 



from which dX^/dcf) can be eliminated to yield the result. Therefore, 



K{4>\d) 



{<t>-my 



lis 



2(1 -d) 2 det(#*(A„,A„ 



giving the Gaussian fluctuations of the FD. The variance as a function of the density 
then reads 



V„U\A\ (l-rf) 2 det(H*(A n (d),0) (1-d) 2 , x y 

4.3 Special case: the TASEP 



(29) 



Let us apply this formula first to the simple M/M/l file corresponding to one single 
speed level (fj, a = fib = ft), i-e. a TASEP process. The rate of arrival A is set by 
convenience to A = dfi. The cumulant-generating function then reads 



h(s,t) = log(l - d) + lof 



l + de- s (e- tv 
1 - de~ s 



1) 
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Therefore we have 



det(#*(0,0)) 



d 



dfx 



(1 - d) 2 

dfi d(l — d)fj, 2 



d 3 a 



1 - d r 



while E[cj>\d] = fid(l - d). Equation (|29f yields 
Var val (0|d) 



1 — ^ j2/-, n2 2 

a (1 — a) (j, . 



L 



(30) 



The direct result, in the grand canonical ensemble, with expectation constraint on the 
density is given by 



Var G c(0|d) = , - ht 



1 - d 



d(l - d). 



L " L 

Actually, the TASEP computation on the ring geometry can be done directly as 

1 1 N+L 

— Var(0|d) = Var(— — ]T Ti f i+ i 



N + L . 

Z — 1 

(L + N)(L + N-3) 



{L + NY 



+ 



1 



L + iV 



E(Tif, 



i+i; 



E(riTi + iTjrj + i) 



(L + iV) 2 



with Ti = 1 — Tt G {0, 1}, tjv+l+1 = T\ and V(i, j) s.t. |i — j| > 2. The expressions 
above are computed as to give 



NL 



and 



^TiTi+iTjTj+i 



(L + iV)(L + iV-l) 

N(N -1)L(L- 1) 



= ' /(| -' /,(l +L^v) + °(Tiv) 



(L + AT)(L + AT - 1)(L + JV - 2)(L + TV - 3) 

d2 ( 1 - d ) 2 ( 1 + ^-^ 1 - d ) + (LTiv)' 



which finally yields ( 30 1 . 



5 Queues with two-state service rates 

Among the possible mappings discussed in Section [2.3| taking empty sites as queues 
should give more information on the jam structure, i.e. presumably on the long range 
correlations of the model associated to cluster formations, than the mapping of type 
(i). Unfortunately, the mapping of type (ii) is not exact in this case. Nevertheless, 
let us try to define a queueing process able to capture these correlations at least 
qualitatively. 

The first step that we take in this direction is to understand the dynamics and 
steady state regime of a single cluster of vehicles of size n t , assuming that vehicles of 
type A [resp. B\ join the queue with rate A a [resp. A&], while they leaves the queue with 

rate /j, a [resp. fib]- Let A = f \ a + At represent the intensity of the incoming process. In 
the bulk of the queue, accordingly to the transitions rules defined in Section [2. 1| only 
deceleration is possible, with rate 5. 
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5.1 Speed profile inside a single cluster 

To study the speed profile, i.e. the speed labels of particle depending on their posi- 
tion inside a cluster, let us encode a given cluster configuration as a binary sequence 
{(Ao, . . . , A n -i}, with Ai = 1 — Bi, £ {0, 1}, corresponding to either a fast or a slow 
vehicle, n represents the size of the queue, the front vehicle has index n— 1 while index 
corresponds to the last entered one. Let us remark first that the front end of the 
queue is a moving interface, which next n moves are completely conditioned by the 
present state of the sequence. Let us forget about this interface for the moment and 
consider the distribution of the infinite sequence {(Ai, A2, A3, . . .}. The joint law has 
not a product form, because the whole sequence is shifted by one unit (Ai — > A i+ i) 
at instants of arrival, which causes dependencies between neighbors. Nevertheless, the 
structure of dependency is quite simple, for example single site marginals are deter- 
mined in close form. Let Pi(t) = Pt(Ai — 1), wc have 

^ = X(j)i-i - Pi) - Spi, Vi > 1 

^=A„(l-pi)-(A fc + *)pi, 

with A = f A a + Xb. The generating function 

00 

Mz) = ^Pi+l{t)z\ 

d(t> t (z) _ 



satisfies 



with solution 



(Az -(A + + A„, 



4>t(z) 



A + S - Xz 



Correlations between nearest neighbors, qi(t) ^= Pt{Ai = 1, A i+ \ — 1) satisfy a similar 
closed form equation: 

dq 



= A(gi_i - qt) - 2Sqi, Vi > 1 



dq 

with 



^ =X a pi(t)-(X + 2S)q ly 



which generating function 



satisfies 

= (Xz - (X + 26))M*) + (*)> 

having for solution: 

A 2 

t/,t(z) = (X + 5)(X + 2S- Xz) 

+ (*•« - (A + *)(A + 2*-As) + V^W* Mm - e ^ t)y 



5 Queues with two-state service rates 



18 




Fig. 6: (a) State flow diagram of a queue with stochastic service rate, (b) Details 
of the transition rates. 



At steady state we see that 

X a I A \ 1 oo Aif / A 



Pi A 



(a + *)' and qr = x(xTT)(xT2sy 

indicating a significant level of correlations when 5/X is O(l). 



5.2 Joined effective process (n(t) , /j,(t)) 

Since the front end interface of the cluster has no causal effect on the rest of the queue, 
except for the front vehicle which may accelerate with rate 7, we can consider the 
dynamics of the sequence independently of the motion of the front interface. Making 
the additional assumption of independence of the local speed labels in the bulk is 
rather crude regarding the results of the previous section, but nevertheless we look 
for a qualitative comparison. To this end let us write a master equation of the joint 
process (n(t) , fJ.(t)) , i.e. the equation governing the evolution of the joint probability 
Pt{n, t) — P(n(t) = n, [i = [i a T + /t/&f), the joint probability that the queue has n 
clients and its front car is of type i (t = 1) or B (f 1 - r = 1). Given 

dot X a n dot A 

Pn = -T-r with r = . 
A A + 

and p n {r) = p n r + p n f, using p n to denote 1 — p n , the master equation then reads 
dPt ^' T) = X(P t {n - 1,t) - Pt(n,T)) + (ti a Pt(n + 1, 1) + fj, b P t {n + 1, 0))p„(r) 
- (fi, a T + fi b f)P t {n,T) + i(t - f)P t {n,0), n> 2 

dPt ^' r) = (X a r + X b f)P t (0) - AP t (l, r) + (ii a Pt% 1) + li b P t (2, 0)) P i(t) 
- (fi a T + Hbf)P t (l, r) + 7 (t - f)P t (l, 0), 

and 

=-AP t (0)+/ io P t (l,l)+/x 6 P t (l,0). 
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It is a special case of a queueing process with a 2-level dynamically coupled stochastic 
service rate, as defined in Section [3] which state-graph is sketch on Figure [6] In the 
stationary regime, we denote 

7r^ = f P(n(t) =n,n = Ha) and 7r* = f P(n(t) = n,H = Hb) 
and 7T = P(n(t) = 0). 

Proposition 5.1. For n > 1 i/ie following recursion holds 



"n+l. 



HaHb + Ha(l + Aj5n+l) + X^bPn+l 

while for n = 1 we have 



7 + + Aj5n+1 7 + Apn+i 

\pn + l Ha + Xp n + 1 



(31) 



XaHb + A7 + A 2 pi 



Ma/if) + ^lMa + ^PlfJ-b + 7Ma 



6 
""I 



XbHa + A pi 



HaHb + XpiHa + XpiHb + "(Ha 



K0, 



7T(). 



Proof. See Appendix [B] ■ 

Note that for large n, because of a finite acceleration rate 7, the effective service 
rate h°° nas a limit slightly above /if,. It is obtained as 



H°° = Hb + r— — (Ha - Hb) 
l + T] 



where 77 is the limit ratio 



77= lim ~ = - 1 -(v / A + 7- \ + Hb-Ha) with A =' (A - 7 + /i a - /i b ) 2 + 4A7 

OO 



n— ^00 7T^ 2A 

obtained from (131 



is automatically fulfilled by virtue of the partial balance equation 1 35). Consider the 
generating functions 

OO 

g a ,b{z) = ^2 TTn b z n and g(z) = f tt + g a (z) + gb(z). 
n=l 

Proposition 5.2. (i) g(z) satisfies the functional equation of the type 

ug(rz) = — (z — z + )(z — z~)g(z) — vz + w (32) 



where z are given by 



and 



= ^{Ha + Hb + X + J±^/A), 



\ 2 U = \ a (Ha — Hb) X 2 V = (XaHa + ^bHb)^0 X 2 W = (/i a Hb + X a Ha + XbHb + 1 Ha )tT0 . 

fii) T/ie solution reads: 



w — vr z 



Wl =0 {zr k - z+){zr k - z )' 



(33) 
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Proof. See Appendix |B| for details. 



The reduced generating functions g a ,b are then obtained using (35 1 to yield 

i \ Mb . A# Mb / \ , ,s Ma - flo , > 
9a(z) = 7T H 5(2), and g^) = 7r H 

Ma — Mb Ma — Mb Mb ~ Ma Mb — Ma 

(34) 

From these it is then straightforward to obtain the 7r^' b upon using Cauchy integrals, 
as sums of geometric laws. 

5.3 Mean field estimates 

If we consider then the closed tandem formed out of these effective queuing processes, 
the first observation is that the sufficient conditions given in section |3.2| may not 
be fulfilled. As such the large deviation estimation of the FD is going to be only 
approximate, a kind of mean field approximation, for a second reason, in addition to 
the simplifying assumption made on the internal cluster structure. 
From ( 34 1 , we write down the cumulant generating function 



h(s,t) = log 7ro(H ) + - g(s) . 

V Ma — Mb Ma ~ Mb 

Therefore we get for the Hessian 
det(#(s*,0)) = 

gg"-g' 2 s(Ae-* ff + (1 - A e - 8 V) 

g(Xe~ s g + (I - Ae~ s )g') g(fj, a ^b^o + (A(/x a + Mb)e~ s - A 2 - ^ a ^b)g) 
where s* is the point which satisfies 



from this we can compute the small (Gaussian) fluctuations of the FD. A comparison 
with the corresponding ABTASEP is displayed on Figure [7] 

In addition to the fundamental diagram, under the canonical ensemble constraint, 
it is also interesting to determine the single queue distribution. This can be obtained 
from the partition function [16] in the large deviation framework as 



pcE(n,n) = tv (n,n) 



Z L (N) 

~ tt x (n, n) exp\L(h(s(d - x),0) - h(s(d),0) d ~~ - s(d-x) + — —s(d))], 

I 1 — d — x 1 — d ' l 

with x = nj (N + L) and the density constraint 

£«*■»- A. 

satisfied by s(d). A comparison of this queueing formulation with the ABTASEP is 
given on Figure [8] In this case we impose an additional self-consistency condition on 
the parameter A a , 

oo 

A a = 7I"nMa, 
n=l 

which otherwise would be free. The correspondence between the cluster size distribu- 
tion observed on the ABTASEP with the single queue distribution obtained from the 
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(b) 

Fig. 7: Comparison between the fundamental diagram simulated with exclusion 
processes and the one obtained from the corresponding (approximate) 
queuing process. The size is L = 1000, 5 is set to zero to have a TASEP 
(a), while the parameters are set to \i a = 10 x = 10 x 7 = 100 x 5 on 
panel (b). 
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(b) 

Fig. 8: Comparison of particle cluster vs queue's size distribution for set of 
parameters /i Q = 10/i;, = 100, 7 = 10 and 5 = 1 for various sizes L with 
fixed density d — 0.35 (a) and various densities with fixed size L = 5000 
(b), the number of queues being (1 — d)L. 
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generalized queueing process is rather accurate. In particular in both cases, a bump is 
observed in the distributions at the same location, for small size systems. It indicates 
that condensation is observed as a finite size phenomena. In the thermodynamic limit 
macroscopic jams are absent. In this respect it is different from the type of condensa- 
tion analyzed in [16], which are obtained under some conditions on the service rate, 
as a large deviation principle but with different scaling (speed in the large deviation 
terminology) then L. Concerning the FD a larger discrepancy is observed between the 
ABTASEP process and its corresponding effective queueing process. The reason for 
this, which is not visible on the cluster distributions, is that the ratio of fast over slow 
particles for given size of cluster size does not coincide. We don't know however which 
approximation between 

• the simplified assumption on internal cluster structure 

• the neglect of correlations between queues 

in the effective model is responsible for that. These results indicate anyway that the 
steady state of the ABTASEP might be well accounted for by a little bit more refined 
joint cluster measure. 

6 Conclusion and Perspectives 

In this paper, motivated by questions raised in the context of traffic modelling, we 
have proposed a simple extension in the definition of the TASEP model to take into 
account acceleration and braking thereby offering the possibility to study the effect of 
asymmetry between the two mechanisms within a simple model. With two different 
speed levels a rich dynamics is already present. An effective mapping on a generalized 
zero-range process allows to interpret at least qualitatively the spontaneous jamming 
phenomena occurring in this model for some choice of parameters. In addition we 
develop a large deviation formalism to study the fundamental diagram associated to 
this family of queuing processes. Still the discrepancy between numerical simulation 
and estimation seen on Figures [8] and [7J may originate from the various approximations 
we make, concerning the detail structure of the jams and the correlations between 
queues which are neglected when using our large deviation estimation of the FD. This 
formalism could be possibly adapted to the situation where the joint measure of queues 
has not a product form. 

In the family of models that we considered we found that condensation phenom- 
ena are associated to finite size effect on the ring geometry, and at least numerically 
the presence of a large macroscopic jam at large scale seems doomed to decay ex- 
ponentially with no bump in the probability distribution, and this is consistent with 
the approximate mapping on zero-range process that we propose. Nevertheless the 
situations could change when the number of speed levels is increased, as shown in our 
previous work |15) on the subject, and we suspect that synchronized flow can take 
place in this limit of large number of speed level. 

The situation where vehicles may enter or leave the system at very low rates so that 
the density may change adiabatically and allow to observed hysteresis effect deserves 
also more studies. In particular it would be interesting to find a relevant measure 
associated to trajectories in the FD, based for example on Brownian windings models 
proposed in [22) . to quantify the hysteresis level in various regions of the FD. 

Finally, the models considered here are limited to single lane traffic and could be 
easily generalized to multi-lane, using coupled exclusion processes like in 23 . 
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A Complement to Section 4.2 constrained partition function 
and the dual Hessian 

In this appendix we explain the role played by the dual Hessian in the Large devia- 
tion expression of the constrained partition functions. Assume that we have n linear 
constraints, written as 

Cy = V, 

where C is n x d matrix, each line C, corresponds to constraints i and V is a n- 
dimensional vector and that we want to estimate the constrained partition function 
Zl[V]. To obtain the small fluctuations we approximate first T in (21 1 at second order 
with respect to some reference point y* 

Hv) = Hvl + (y - y') T ■ vHy") + (y- y*) T v 2 Hv'Xv - y*) + °(l|y - y* II 2 ) , 

so that y considered as a d- dimensional vector is integrated over R d . y* is chosen such 
that 

y* (A) = argmin \T[y\ - A T Cy , 

y L J 

matches the constraints with the proper Lagrange multipliers A. Since in absence of 
constraints Zl[V], as it stands in ( |20[ ) is normalized to 1, so at this order of approxi- 
mation the partition function reads 

Z L [V]=L d/2 detH f d d y e -L(T(y*) + (y-y*) T -vF(y*)+^y-y*) T H(y-y>)) 
V (27I~) d J 

n d 
i=l 3=1 

with 

h = v 2 Hvl, 

the Hessian taken at y = y* . The dual Hessian is associated to the dual free energy 

T*[A]=T[y*(A)]-A T Cy. 
Using relations associated to the stationarity of y* , it reads 

H* = C T H^C. 



Due to the specific form (22 1 of _F, the Hessian turns out to be the covariance matrix 
between the various quantities Vt, as given in (281. A convenient way to express the 
constraints is to write 



so that a Gaussian integration over y can be performed and since 

\\V-Cy\\ 2 = \\C( y -y*\\ 2 , 

we simply get 



n^oll (2na) n det H a 
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with 

La 

Let P be the orthonormal projection on the subspace spanned by the set of constraint 
vectors Cu an P = f 1 — P, such that 

C = CP and P C T = C T . 

We have 

det(fla) = det(#)det(l + ^-^H' 1 ) 

Lot ' 

= detfi?) det(P(l + C ^-H- 1 )P + P C - C ^H~ X P + PP) 
La La ' 

= det(H)det(l + ^^H' 1 ) 
p v La ' 

where detp corresponds to the block determinant associated to the subspace of con- 
straints, and 

H* = C T H^C 



is the dual Hessian defined in (28 1. Concerning V a , since y* is a stationary point at 



least w.r.t. fluctuations orthogonal to the constraints, so we have 

PV a = 0. 

Since 

v La ' 
= H' 1 (PLaH(C T C)~ 1 P — PP + PP) + o(a), 

as a result 

lim V^H' 1 = 0. 
So finally the Gaussian estimate of the partition function reads 

rn/2 

Z L [V] = , L e~W> 
x /(27r)"det( J ff*) 

which yield in particular the expressions ( |26|27[ ) up to a factor L n caused by a different 
convention in the constraint definition. 

B Complement to Section [5j steady state of the queueing 
process 

From the master equations, taken at steady-state, we get the recurrence for n > 2 

\ / a a\ . { a . b \ a , b r\ 

A[Tr n -l — K n ) + {fl a TT n + l + ^b^n+ljPn — ^a^n + l^n = 0, 

\ f b b \ , / a i b \ - b b n 

A(7T n _l ~ 7T n j + (fJ, a TV n + 1 + ^ 6 7T n + 1 jp n - fl b TY n - JTT n = 0. 
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For n = 1, it reads 

X a TV - XlT 1 + [fl a ^2 + f-b^jPl ~ Ma^l + J^l = 0, 
Ai,7T — Xni + {fla^l + flb^Pl — fi-b^X — 1^1 = 0. 

The sum of the two equation gives that the quantity Xir n — {u a 'Kn+\ + f^b^n+i) i s a 
constant independent of n which has to vanish, leading to the partial balance relation 



Using this gives the relations, 

(fi a + Ap n+ i)< +1 - (7 + Ap„+i)7r£, +1 = A< 

-Ap„+l7T, a i + 1 + (7+ fl b + \p n + 1 )% b n + 1 = \-K h n . 



(35) 

(36) 
(37) 



The recurrence is then inverted to yield (31 1. 

Recurrence ( 36|37 1 are then conveniently used to obtain the following relations 
among the generating functions, 

X a g a (rz) + \ a gb{rz) = (pL a + A - Xz)g a {z) — 7ffb(z) — X a woZ 

X a ga(rz) + X a gb{rz) = \g a (z) — (7 + fib — \z)g b {z) + XfKoz. 

It follows that 

. . Xz \z — fx a . . 

9b(Z) = r-7T + r JaW, 

fib — Xz fib — Xz 



also related to ( 35 I , which rewrites 
Xz — fib 



9a(z) 



fib 
fla — fib 



-n + 



fla — fib 



and 



9b{z) 



fla , Xz fl a / \ 

7T H g(z). 



fib — fla 



fib — fla 



So finally g(z) satisfies the functional equation given in proposition 5.2 The solution 
is constructed as follows. Consider the infinite product 



G(z)^H(l 



r)(l 



We have 



G{rz) 



so a solution of the form 



(z — z ){z — z+) 
C(z 



+,G{z) 



(38) 



G(z) 



has to verify 

, U „, s W — VZ 

Cz + -+^C(rz) = — 

z + z (z — z ){z~z + ) 

Equivalently, for any n > this rewrites 



._^_)»C(r-,)-(-^ r )- +1 C(r^ 1 ,) = (- 



G(z) 



w — vzr 



z + z (zr n — z )(zr n — z + ) 



= {-u) r 



w — vzr 



nLo( zrfc - z )( zrk - z+ ) 



G{r n z) 



G(z) 



after using relation ( 38 1 . Taking the sum leads to the solution ( 33 1 , owing to 

lim (- -^—) n C(r n z) = 0, 
which can be checked afterwards. 
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